Bardeen variables and hidden gauge symmetries 
in linearized massive gravity 
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We give a detailed discussion of the use of the (3 + 1) decomposition and of Bardeen's variables in 
massive gravity linearized over a Minkowski as well as over a de Sitter background. In Minkowski 
space the Bardeen "potential" <!>, that in the massless case is a non-radiative degree of freedom, 
becomes radiative and describes the helicity-0 component of the massive graviton. Its dynamics is 
governed by a simple Klein-Gordon action, supplemented by a term (D$) 2 if we do not make the 
Fierz-Pauli tuning of the mass term. In de Sitter the identification of the variable that describes 
the radiative degree of freedom in the scalar sector is more subtle, and even involves expressions 
non-local in time. The use of this new variable provides a simple and transparent derivation of the 
Higuchi bound and of the disappearance of the scalar degree of freedom at a special value of m 2 g /H 2 . 

The use of this formalism also allows us to uncover the existence of a hidden gauge symmetry of 
the massive theory, that becomes manifest only once the non-dynamical components of the metric 
are integrated out, and that is present both in Minkowski and in de Sitter. 



I. INTRODUCTION 

The problem of formulating a consistent theory of grav- 
ity with a massive graviton has a long history. Already 
in 1939 Fierz and Pauli [1] showed that, when the the- 
ory is linearized over Minkowski space, a specific form 
for the mass term is required to avoid the appearance 
of a sixth ghost-like degree of freedom. In the 1970s, 
in another classic paper, Boulware and Deser [2] showed 
that even with the Fierz-Pauli (FP) mass term the ghost 
reappears when the non-linearities of the gravitational 
field are taken into account. Recent years have witnessed 
a flurry of activity on the problem, stimulated by the 
possible relevance of infrared modifications of gravity for 
understanding the origin of dark energy, and by theo- 
retical breakthroughs [3-5] that also have an intrinsic 
field-theoretical interest. 

A recurrent theme in the study of both gauge theo- 
ries and general relativity is that a good choice of vari- 
ables can significantly enlighten the physics, and different 
choices are appropriate for studying different aspects of 
the theory. Indeed, the recent breakthroughs in massive 
gravity can be partly traced to a clever way of isolating 
the dynamics of the helicity-0 mode through a general- 
ization of the Stiickclbcrg trick to general relativity, pro- 
posed in [3] (and further studied in various contexts in 
[6-10]. See [11] for a recent review). This lead recently to 
the construction of a consistent ghost-free theory of mas- 
sive gravity to all orders in the decoupling limit and up to 
quartic order in the non-linearities away from the decou- 
pling limit [4, 5, 12]. The absence of ghosts in this model 
was proved in full generality using the ADM formalism 
in [13]. Ghost-free actions with a general reference met- 
ric that admits massive spin-2 fluctuations around non- 
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flat backgrounds were first considered in [14], and were 
proven to be free of the Boulware-Deser ghost in [15-17]. 
Ghost-free dynamics for the reference metric is presented 
in [18]. Further recent work on ghost-free bimetric theo- 
ries includes [19-22]. 

In this paper we provide a detailed discussion of the 
use of the (3+1) decomposition of the metric and of the 
gauge-invariant Bardeen variables in linearized massive 
gravity. This formalism is a standard tool of cosmological 
perturbation theory [23]. In massive gravity linearized 
over Minkowski space it was first introduced in [24] , and 
has been applied to massive gravity in a number of recent 
papers, sec e.g. [25-27]. Elaborating on these results, we 
will see that these variables can be quite useful for eluci- 
dating various aspects of the massive theory. In particu- 
lar we will see that in Minkowski space, as the mass term 
approaches the FP form, after elimination of the non- 
dynamical components of the metric the dynamics in the 
scalar sector nicely collapses to a massive Klein-Gordon 
action for the Bardeen variable $, which therefore de- 
scribes the hclicity-0 component of the massive graviton 
(see also [24, 26]). It is quite interesting to see that $, 
which in the massless case describes a non-radiative de- 
gree of freedom (and is in fact usually called the Bardeen 
"potential" ) becomes a radiative degree of freedom when 
the FP mass term is switched on. Furthermore, for a 
generic mass term the action in the scalar sector can be 
reduced to a higher-derivative theory for $, in which only 
Lorentz-covariant structures such as (D$) 2 appear. More 
generally, we will see that the use of the Bardeen variable 
$ can be a convenient way of isolating the dynamics of 
the helicity-0 mode, complementary to the by now stan- 
dard Stiickclbcrg formalism. 

An intriguing consequence of the fact that the helicity- 
mode of the massive graviton is described by $ is 
that a residual gauge symmetry appears, consisting of 
the transformations — > — {d^ v + d^^) with £ M 
parametrized by two scalar functions A(x) and C(x) as 
£o = A and & = diC . We will see that this symme- 
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try only appears when one eliminates the non-dynamical 
components of the metric, remaining just with the five in- 
dependent fields that describe the physical components 
of a massive graviton (plus the extra ghost-like scalar if 
we are away from the FP point). 



We will then turn to massive gravity linearized over 
a de Sitter background. We will see that in this case 
the identification of the radiative degree of freedom that 
describes the helicity-0 mode of the massive graviton is 
quite subtle, and we will show that it even involves in- 
tegrals over time of some metric components. In terms 
of this variable, after elimination of the non-dynamical 
degrees of freedom, the dynamics in the scalar sector of 
linearized massive gravity in de Sitter (with FP mass 
term) collapses again to a simple Klein-Gordon action. 
The overall coefficient of this action nicely displays the 
existence of the Higuchi bound as well as of a point of 
enhanced symmetry corresponding to the so-called "par- 
tially massless" theory, where the scalar degree of free- 
dom disappears. At the same time this dynamical vari- 
able is gauge-invariant, which shows that even in de Sit- 
ter space, after elimination of the non-dynamical fields, 
gauge transformations with £o = A and = diC are a 
symmetry of the massive theory. 



The paper is organized as follows. In sect. II we recall 
the (3+1) decomposition of the metric in flat space-time 
and we illustrate how it can be used to separate the met- 
ric into pure gauge, radiative and non-radiative degrees 
of freedom. We also discuss the behavior under Lorentz 
transformations of the variables entering the (3 + 1) de- 
composition and of the gauge-invariant Bardeen's vari- 
ables. Most of the material in this section is known in 
the literature (except for the part on Lorentz transforma- 
tions of the variables entering the (3+1) decomposition), 
but we find useful to present in a systematic way various 
results that will be needed in the rest of the paper. In 
sect. Ill we use these variables to study massive gravity, 
linearized over Minkowski space, for a generic quadratic 
mass term. We explicitly identify the ghost degree of free- 
dom and we find that, at the FP point, the scalar sector 
is described by the Bardeen variable We will then 
see that, outside the FP point, the scalar sector can be 
reduced to a simple higher-derivative theory for $. The 
hidden gauge symmetry that emerges from this analysis 
is discussed in sect. IV. In sect. V we discuss this symme- 
try from the point of view of the Stuckelberg formalism. 
In sect. VI we compare our results with a similar analy- 
sis performed in massless and massive electrodynamics. 
In sect. VII we discuss massive gravity linearized over 
a de Sitter background. Sect. VIII contains our conclu- 
sions and a summary of the main results. Some technical 
material is relegated in appendixes. We use the signa- 
ture r\^ v = (— , +, +, +), units h = c = 1, and we define 
k = (327TG) 1 / 2 . 



II. DECOMPOSITION OF THE METRIC INTO 
PURE GAUGE, RADIATIVE AND 
NON-RADIATIVE DEGREES OF FREEDOM 

A. (3+1) decomposition of the metric 

We begin with a review of the (3+1) decomposition in 
the massless case. Expanding the metric g^ v over flat 
space, g^ v = rj^ v + h^, the metric perturbation can be 
decomposed as 

h 00 = 2ip, (1) 

h 0i = ft + 5,7 (2) 

hij = -20% + (didj - ^SijV 2 ) A 

+\{d i e j + d j e i ) + h^ , (3) 

(where V 2 is the flat-space Laplacian) with the con- 
straints 

a,ft = 0, 8^ = 0, Wig T = 0, S ij hj^ = 0, (4) 

and the boundary conditions 

7^0, A^O, V 2 A^0, e 4 ^0, (5) 

at spatial infinity. As we will review below, these bound- 
ary conditions ensure the uniqueness of the decompo- 
sition [28]. The tensor hj? is a transverse traceless 
(TT) 3x3 symmetric tensor, so it carries two degrees 
of freedom. The variables (3 l and e 1 are transverse vector 
fields, so they carry two degrees of freedom each, and 
describe vector perturbations of the background. Fi- 
nally, we have four fields ?/>,</>, 7, A that are scalars un- 
der spatial rotations and describe the scalar perturba- 
tions of the background, for a total of 10 degrees of free- 
dom. This parametrization decomposes the metric per- 
turbations into irreducible representations of translations 
and spatial rotations and corresponds to a decomposition 
into eigenfunctions of the Laplacian, i.e. to harmonic 
analysis. For this reason it is also called the harmonic 
decomposition, and we will also refer to the variables 
{ip, </>, 7, A, (3 1 , e i , hJ } T } as harmonic variables. 

The linearized massless theory is invariant under the 
gauge transformations 

(d^ v + d„£ M ) , (6) 

which corresponds to linearized diffcomorphisms. To un- 
derstand the properties of the harmonic variables under 
this gauge transformation it is useful to write the gauge 
functions £ M in the form 

£o=A, Z i = B i + d i C, (7) 

where Bi is a transverse vector, d{B l = 0. In terms of 
these variables eq. (6) reads [28] 

i)^i>-A, 0^0+Iv 2 C, (8) 

j^j-A-C, A^A-2C, (9) 
A -+Pi-Bi, e^e,- 2B t , (10) 
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while hj^- is gauge invariant. As dictated by symmetry, 
the transformation of the scalars ip, </>, 7 and A depends 
only on the scalar functions A and C, while the trans- 
formation of the transverse vector fields Pi and ti only 
depends on the transverse vector field Bi. The fact that 
hj^~ is gauge invariant is a consequence of the fact that, 
from the point of view of spatial rotations, £ M decom- 
poses into a spin-0 and a spin-1 part, while a traceless 
symmetric tensor such as hJ T is a spin- 2 operator. 

Using the above variables one can form the following 
gauge-invariant scalar combinations 



that A itself vanishes at infinity allows one to invert once 
more the Laplacian in eq. (16) and obtain A. From these 
expressions for and V 2 A we find that the Bardeen vari- 
able $ can be written as 



4$ = h\- v^afr&hi 



(17) 



Thus, even $ is a non-local function of /i M „, and the same 
holds for <5, which involves A and therefore a double in- 
version of the Laplacian. In the vector sector the inver- 
sion of the harmonic decomposition gives 



-t-\v 2 x. 



(11) 

(12) 



whose generalization to perturbations of FRW space-time 
gives the standard Bardeen variables. We will still use 
this nomenclature in the flat-space case. Similarly we can 
form a gauge-invariant transverse vector 



5, = A 



1 



(13) 



Thus, at the level of linearized theory, we have six gauge- 
invariant quantities: the two components of the trans- 
verse traceless tensor perturbations hj? , the two compo- 
nents of the vector perturbation 3j subject to the con- 
dition diS 1 = 0, and the two scalar perturbations $ and 
So, the four gauge functions £ M allow us to eliminate 
four pure-gauge degrees of freedom from the ten compo- 
nents of h^y , remaining with six gauge- invariant degrees 
of freedom. 

It is important to appreciate that, in a generic gauge, 
the harmonic variables are in general non-local functions 
of the metric [28]. This can be seen inverting eqs. (1)- 
(3), as follows. The variable ij> is simply given by eq. (1), 
while </> is obtained taking the contraction of eq. (3) with 
6 ij , 



1 



'00 : 



(14) 



Thus, these quantities are local functions of the metric. 
All the other variables, in contrast, have a non-local de- 
pendence on hoi or on hij. To extract 7 we take the 
divergence of eq. (2) and we invert the Laplacian (which, 
with the boundary condition that 7 vanishes at infinity, 
is a well-defined operation). This gives 



7 = V~\d l h 0l ) . 



(15) 



To extract V 2 A we apply the operator d l d^ to eq. (3) 
and we get 



ft = hoi-V-^di&hoj), (18) 
U = 2V- 2 [&h lJ -d l V- 2 (d k d l h kl )] , (19) 



while hj? can be obtained from eq. (3), using the above 
expressions for 0, A and e,-, and involves both the operator 
V~ 2 , and the operator V -2 V -2 . Observe however that 
in the massless theory, when T^ v = 0, we can use the 
gauge invariance to fix the TT gauge, where h\ = 0, 
hop = and 9*/iy = 0. In this gauge the scalar and 
vector variables vanish, as well as the non-local terms in 
hJ T , and = hJ T . 



B. Radiative and non-radiative degrees of freedom 



1. Action and equations of motion for the gauge-invariant 
variables 



To study the dynamics of the harmonic degrees of free- 
dom we consider the linearization of the Einstein action, 
expanding g^ v = rj^ u + h^ v . The quadratic part of the 
Einstein-Hubert action and the interaction term with an 
external conserved energy-momentum tensor are given 
by 



S2 + Si n t 



d A x 



2/s 2 ^ 



1 



(20) 

where n = (32ttG) 1 ^ 2 and £^ V 'P" is the Lichnerowicz op- 
erator, defined as 1 



+ Wd^d" + ri" v d p d a ) , 



(21) 



and □ = rj^dpdv is the flat-space d'Alembertian. It is 
useful to perform the harmonic decomposition also in the 



V 2 A 



V~ 2 (d l d J h l:j ) . 



(16) 



where we used the boundary condition that V A van- 
ishes at infinity to invert V 2 (V 2 A). Requiring further 



1 Beware that, in the literature, different conventions are used for 
the overall sign of the Lichnerowicz operator. With our conven- 
tion £»»<P u h P a = +Uh» v - .... 
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energy-momentum tensor, writing [28] 



The result is 



Too 
Toi 

Ti- 



ft, 

Si + d t s , 



PS, 



didj - -5 l3 V z ) a 



-{d l a + djcr t ) + Ui 



(22) 
(23) 



(24) 



5 2 + S h 



d 4 x 



12$ 2 + 4<9 l $<9 ? $ - 8c\$<9 1 * 



J d 4 x 3$P + ^ P - 5,5* + ^hj?a ij 



(31) 



where 



d t a l = o , d t s l = o , aVi 



0. 



0. (25) 



The isotropic part of is P#y, where P is the pressure. 
The remaining terms in Ty define the anisotropic stress 
tensor and depend on a scalar a, a transverse vector Oi 
and a TT tensor cry. As with the harmonic decomposi- 
tion of the metric, the uniqueness of the decomposition 
is assured if we assume that 5, a, V 2 cr and <t ? vanish at 
spatial infinity. 

The quantities that appear in the paramctrization of 
7)t„ are not all independent, since they are related by 
energy-momentum conservation d^T 1 *" = 0. Imposing 
d^ = gives 



V 2 5 



while d^T^ 1 = gives 

i V V - S* + d i f-S + P + | V 2 <t 



(26) 



(27) 



As expected, the action depends only on the gauge- 
invariant combinations Si and hJ 3 T . We can now 
derive the equations of motion. Using {4>,ip, A, 7} as in- 
dependent variables, the variation with respect to <j) gives 



24$ - 8V 2 $ + 8V 2 * = -3k 2 P . 



(32) 



while the variation with respect to ^ gives the Poisson 
equation 



8V 2 $ = -n 2 p. 



(33) 



The variations with respect to A and 7 give a combi- 
nation of derivatives of these equations. Observe that 
eqs. (32) and (33) can also be obtained taking directly 
the variation of the action (31) with respect to and *, 
respectively. 

The variations with respect to /3i (or, equivalently, with 
respect to 3*) and to hJ T give, respectively 



2V 2 5. ; , = - K 2 S t 



2Uh. 



TT 



-K Oi 



(34) 
(35) 



Since o l and S % are transverse, see eq. (25), applying di 
to this equation we get the condition 



V 2 [-5 + P + (2/3) V 2 a] = 0. 



(28) 



We impose the boundary conditions that the energy- 
momentum tensor vanishes at infinity (or just that 5, 
P and V 2 cr vanish at infinity). Using the fact that a 
Poisson equation V 2 / = with the boundary condition 
/ = at infinity only has the solution /(x) =0, we get 



-5 + P + -V 2 ct = 0. 



Inserting this into eq. (27) we then get 



vV 



25' 



(29) 



(30) 



Observe that the vector equation (27) separate into an 
equation for the transverse vector part and one for the 
scalars that parametrizes the longitudinal vector part. 
Thus, overall energy-momentum conservation gives two 
scalar conditions, (26) and (29), and one condition be- 
tween transverse vectors, eq. (30). 

We can now write the linearized action using the (3 + 1) 
decomposition of the metric perturbations, eqs. (l)-(3). 



while the equation obtained performing the variation 
with respect to Ci is implied by eq. (34). Observe that 
* enters linearly in the action, so it is a Lagrange mul- 
tiplier. Integrating by parts the term — 8di&d ll $> in the 
action, the part of the Lagrangian that depends on * is 
[(8/k 2 )V 2 $ + p]* and the variation with respect to 
enforces the constraint (33). 

Plugging eq. (33) into eq. (32) we can rewrite the latter 

as 



$ + _V 2 v]/ = —(p + 3P) 



(36) 



The term <3> can be eliminated observing that eq. (33) 
implies that 8V 2 $ = —n 2 p. Using eq. (26), this becomes 
8V 2 $ = -k 2 V 2 5. In flat space it is natural to impose 
the boundary condition that $ and 5 vanish at infinity. 
Using again the fact that a Poisson equation V 2 / = 
with / = at infinity only has the solution /(x) = 0, we 
get 

8$=-k 2 5. (37) 
Therefore eq. (32) can be rewritten as 

V 2 * = -4ttG{ P + 3P - 35) , (38) 
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which can be further simplified using eq. (29) to write 



3P - 3S 
we have 



-2V cr. In conclusion, in Minkowski space 



V 2 $ = 
V 2 * = 
V 2 - = 



-4nGp , 

-4irG(p - 2V 2 ct) , 
— IGirGa.j j . 



(39) 
(40) 
(41) 
(42) 



We see that only the tensor perturbations obey a wave 
equation. The gauge-invariant scalar and vector pertur- 
bations obey a Poisson equation, and therefore represent 
physical but non-radiative degrees of freedom, which are 
fully determined by the matter distribution. We further 
observe that 



V 2 ($ = -8ttGVV 



(43) 



Therefore, if the scalar part of the anisotropic stress ten- 
sor vanishes, we have $ = VP. In the absence of matter 
we have $ = = and Hj =0, and we only remain with 
the two radiative degrees of freedom described by hj? , 
i.e. with the two polarizations of the massless graviton. 



2. Lorentz invariance in harmonic variables 

As with any choice of variables in a theory with gauge 
invariance, the use of the set of harmonic variables 
{ip,(j>, 7, A, /3\e 4 ,/i^ T } has some advantages and some 
drawbacks. The main advantage is that, out of them, we 
can construct quantities which are invariant under lin- 
earized gauge transformations. Furthermore, under spa- 
tial rotations the transformations of the harmonic vari- 
ables are simple: ip, A, 7 and (j) are scalar, P % and e* are 
vectors and hj? is a tensor. In contrast, the behavior of 
these variables under Lorentz boosts is quite complicated 
(and, as we will see below, even non-local). For instance, 
we know that the action (31) is Lorentz invariant, since 
it is just a rewriting of eq. (20) with different variables. 
However, in the form (31) the Lorentz invariance of the 
theory is not at all evident. Furthermore, the scalar, 
vector and tensor sectors are not separately Lorentz in- 
variant, since of course Lorentz boosts mix scalars with 
vectors, while vectors mix with scalars and tensors, and 
tensors mix with vectors. It is therefore interesting to 
study in some detail how the harmonic variables trans- 
form under Lorentz boosts, and how the invariance of the 
the action (31) comes out in this formalism. 

A Lorentz transformation of a tensor can be decom- 
posed in two parts, the spin part mixing the Lorentz 
indices p and the orbital part mixing the field indices 
x^ . It is the spin part that will mix the harmonic vari- 
ables among themselves since we have decomposed the 
Lorentz indices fi — > (0, i). From the point of view of the 
orbital part, these are just ten field representations of 
the Lorentz transformations. This is to make clear that 
the formalism will lose the manifest Lorentz covariancc 



only as far as the mixing of the tensor components is 
concerned. 

So consider an infinitesimal Lorentz transformation 



A>*„x u with A^ 



6Z 



compute the variation Sh^ defined as 



Sh,, 



We then 



(44) 



so that the orbital part does not appear since we are 
only interested in how the components mix. Using the 
expression of the harmonic variables in terms of the met- 
ric found in sect. II A we can obtain their variation un- 
der Lorentz transformation. We restrict to boosts, i.e. 



0, ojq' 1 7^ 0. For tp and <p we get 



Sip 
35cj) 



-tV(A + ^7) 



(45) 
(46) 



Observe that ip + 30 is Lorentz invariant, as it should 
since h = r/^h^ = —2(ip + 30). Furthermore, these 
transformations are local functionals of the metric, since 
(Pi + di"/) is simply h^. The transformations of 7 and 
V 2 A are more complicated and non-local, since their ex- 
pression in terms of the metric involves the inverse Lapla- 
cian. Rather than dealing with the transformation under 
boosts of the inverse Laplacian it is convenient to start 
from V 2 7 = c^/ioi, and observe that, if under a Lorentz 
transformation a quantity / has a variation 5f, then 

S(d IM f)=u lt v d v f + d lt Sf, (47) 

i.e. we treat formally as any other four-vector, satis- 
fying 5(8^) = u^dy. In particular, under boosts 

6(d ) = uj 0i di , 5{di) = ujQido , (48) 

and therefore 

S(V 2 f) = 2oj Q l dJ + V 2 (5f). (49) 

Thus <5(V 2 7 ) = 2w i d i7 + V 2 (<h), while S(d l h oi ) = 
<jJo l hoi + d z 5hoi, and therefore 

V 2 (<S 7 ) = -2to Q l da + u Jo l h 0l + d l Sh 0l . (50) 

This gives 



c>7 = ojq 1 V~ 



0, 



2^-7-20+ -V 2 A 



Pi + ^ V 2 e, 



Using eq. (16) and proceeding similarly, we get 



2<9i7 + -e, 



5{V Z \) = u l -P 



(51) 



(52) 



while in the vector sector, using (18) and (19), we find 



8p % = -u Q 3 V 



5^(2^-20-7) + ^ 



(2iJj - 20 - 7)^ + - (didj - h 3 V 2 ) A 



2 3 



lTT 



(53) 
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Set = 2u j V~ 2 
1 



(<3^-%V 2 )(A- 7 ) 



; tt 



(54) 



For the gauge-invariant combinations we get 

5$ = ^U! %, 



(55) 



<5* 



(56) 



53{ = w J V 



□^T_ (9i H i + a iH .) 



-2(^ -c%V 2 )($ + *) . (57) 



TT 



-CJ V 



'<>J'V I <>j,Y) 



(58) 



Observe that (3 n = (3 l + 5(3 l is transverse with respect 
to the transformed coordinate x n , i.e. d(3 n jdx' % = or, 
equivalently, S(di/3 l ) = (and similarly for ej, Ej and 

Eqs. (55)-(58) show that under boosts the gauge- 
invariant variables transform among themselves, al- 
though their transformation involves the inverse Lapla- 
cian. 2 It is instructive to check explicitly that the action 
(31) is indeed Lorentz invariant. We neglect for simplic- 
ity the interaction term and we split the various terms in 
cq. (31) into the scalar, vector and tensor sectors, 



-^2, scalar 
•^2, vector 

^2, tensor 



-12$ 2 + 4cW$- 



, (59) 

(60) 
(61) 



Under spatial rotations £ 2 , scalar, £ 2 , vector and £ 2 , tensor 
are of course separately invariant. Under boosts, using 



2 Since under rotations <5<t> = Sty = 0, <5Hi = uv'Sj and ShJ^ = 
w i k ^kj 1" w i k ^ki > these variables transform among themselves 
under the full Lorentz group. Of course, one should not be sur- 
prised to find six objects that transform among themselves un- 
der the Lorentz group, despite the fact that symmetric traceless 
tensors form a 9-dimensional irreducible representation of the 
Lorentz group. The point is that these are representations on 
fields, and therefore are infinite-dimensional. This is completely 
analogous to the fact that a transverse vector field fi*(x) satis- 
fying di0* = is a representation of the rotation group, despite 
the fact that it has only two independent (field) components. 
In other words, a transverse vector field is a representation of 
50(3) of dimension (oo) 2 , despite the fact that a vector is an 
irreducible representation of dimension 3. 



eqs. (55)-(58) (and neglecting total derivatives) we get 

^scalar = +WS< V 2 ($ + ¥) , (62) 

^2,vector = -4u> %V 2 ($ + *)-2Lj i E :i Uh% r ,(63) 



5£ 2 ,tensor = +2u i E. j nhJ? , 



(64) 



so the total action is indeed invariant. In contrast 
£2,scaiar, £2, vector and £2, tensor are not separately invari- 
ant (unless we impose the equations of motion for $, 
and tijj 1 - , which in the case = that we are consid- 
ering read V 2 $ = V 2 * = □/i^ T = 0). 



III. MASSIVE GRAVITY IN HARMONIC 
VARIABLES 

We next use these variables to discuss the massive 
theory linearized over Minkowski space. It will be in- 
structive to work with a generic Lorentz-invariant mass 
term, rather than specializing to the FP combination 
from the beginning. It is also convenient to expand 
9nu = V^v + K Vi where n = (32-kG) 1 ^ 2 , i.e. to replace 
h^, v — > Kh^v in the formulas of the previous section, so 
that henceforth has canonical dimensions of mass. 3 
We therefore consider the Lagrangian density 

£2 = \ [V^'^V - + b 2 h 2 )] . (65) 

The Pauli-Fierz point corresponds to bi = 1, b 2 = — 1. 
The equations of motion in the absence of matter are now 



EM h pa = rrPJhhr + b 2 ^ v h) . 



(66) 



Since d^^^' pa h pa ) = 0, they imply the conditions 



duibih^ + b 2 r ) flv h) =0. 



(67) 



A. Elimination of the non-dynamical fields. Scalar 
sector 

We first consider the scalar sector of the theory. We 
begin by writing the action in terms of the variables 
<f), ijj, A, 7 defined in eqs. (l)-(3). We get 



£ 



2,scala 



= -12$ 2 + 4a i $ ( 3 i $-8<3 i *<9 l $ 

h Ui> 2 + 120 2 - 2<9 ? ;7<9 l 7 + ^(V 2 A) 2 



-4&aW> + 30) s 



(68) 



3 Thus, henceforth <j>, ip, <E> and Vt 1 have dimensions of mass, fi is 
dimensionless and A has dimensions of [mass] -1 (while, when we 
perform this decomposition on the dimensionless metric, <j>, ip, $ 
and <E> are dimensionless, 7, = [mass] -1 and A = [mass] -2 ). 
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Since the mass term breaks gauge invariancc, the La- 
grangian now depends on the four fields <j>, ip, A, 7, rather 
than just on the two gauge-invariant combinations $ and 
"J. We now want to eliminate the non-dynamical fields 
and identify the variables which describe the radiative 
degrees of freedom in the scalar sector. A non-dynamical 
variable is integrated out by using its own equation of 
motion if it enters the action quadratically, or by us- 
ing a constraint imposed by a Lagrange multiplier. Such 
an elimination procedure goes through even at the path- 
integral level. The first step is therefore to chose a con- 
venient set of independent variables. 4 



1. Outside the FP point: 61 + 62 7^ 0. 

We find convenient to use <5>,ip, A, 7 as independent 
fields. Observe that the change of variables from <j> to 
$ does not involve time derivatives and is thus legiti- 
mate, see footnote 4. We discuss separately the case 
bi + b 2 7^ and the case b\ + b 2 = (which includes the 
FP point 61 = 1, 62 = — 1 and its sign- reversed), since the 
structure of the non-dynamical equations is different. We 
begin with the case b\ + b 2 7^ 0. 

By inspection of the Lagrangian (68) one immedi- 
ately sees that 7 and ip are non-dynamical [25]. In- 
deed, taking the variation with respect to 7 (and recall- 
ing from eq. (12) that \1/ in eq. (68) contains 7, since 
* = ip - 7 + (1/2) A) gives 



Here we see that the Pauli-Ficrz mass term, for which b\ + 
b 2 = 0, is special. When b% + b 2 7^ we can use eq. (72) 
to eliminate even ip from the Lagrangian, remaining with 
a Lagrangian C 2 {<&, A). As expected, we could eliminate 
two non-physical fields and remain with two degrees of 
freedom in the scalar sector. 

Actually, even A could be integrated out. In fact, tak- 
ing the variation with respect to V 2 A gives 



1 



1 



1 



-m 2 (&i+o 2 )V"A = <t>--m 2 g {b 1 +3b 2 )<S>+-m 2 g b 2 iP . (73) 

Thus, for 61 + 62 7^ we can eliminate V 2 A in favor of $ 
and $ using eq. (73) (with ip expressed in terms of V 2 A 
and $ through eq. (72)). This results in a Lagrangian 
that depends only on However such a Lagrangian in- 
volves higher-derivative terms proportional to ($) 2 , aris- 
ing from the terms (V 2 A) 2 in eq. (68). A theory whose 
equations of motion are fourth-order in the time deriva- 
tives propagates twice the number of degrees of freedom 
than is apparent from its field content. This is due to the 
fact that to evolve the classical equations of motion we 
need to specify twice as much initial conditions [11, 29]. 
Furthermore, one of these two degrees of freedom is a 
ghost [9]. 5 Thus, if we eliminate V 2 A, the corresponding 
(ghost-like) degree of freedom does not disappear, but 
just hides in the fourth-order equations of motion for 

For the moment, rather than eliminating V A, we 
trade it for a new field T defined by 



V 2 (4$-6im 2 7 ) =0. 



(69) 



With the boundary condition that $ and 7 vanish at 
infinity this equation is equivalent to 



4$ - Oim 2 7 = . 



(70) 



m 2 a b 2 



V 2 $ - 12m 2 6 2 $ - 2m 2 6 2 V 2 A) , (74) 



and we use ($, T) as independent variables. The variable 
T is defined so that eq. (72) reads 



In the massless case this reduces to $ = 0, which is in 
fact the same as eq. (37), since we have set = 0. 
When m g 7^ 0, we can rather use it to eliminate 7 from 
the action, 



7 



b\m 2 



(71) 



Thus, in the massive case 7 is a non-dynamical variable 
that can be eliminated algebraically. The variation with 
respect to ip gives another algebraic equation, 

4m 2 (oi + b 2 )iP = 8V 2 $ + 12m 2 o 2 $ + 2m 2 6 2 V 2 A . (72) 



1> = - 



8( 0l + b 2 ) 



r. 



(75) 



Eliminating 7 and ip with the help of eqs. (71) and (75) 
and expressing V 2 A in terms of T and $ we finally get 

Scalar ($,T) = 12$ 2 + 4 ^ ± ^- 8^3% - 12m 2 6i$ 2 

, . . _ 16(6 1 + 02) I ■ _ h(b 1+ 2b 2 ) 

32(h + b 2 ) 9 



b\b 2 m? g 



(76) 



4 Observe that not any choice of variables is legitimate. In par- 
ticular, the initial conditions on the metric hij and on its time 
derivative hij must be in bijection with the initial conditions 
on the new set of variables. For instance this is the case for 
the harmonic variables, but it is not the case for combinations 
of variables involving e.g. A or A, such as This point will 
be important when we consider massive gravity linearized over 
de Sitter, in sect VII. 



° Adding to the mass term higher-order polynomials in /i M „ with 
appropriately chosen coefficients one can have a well-defined 
Cauchy problem with second-order equations and no ghost, see 
[4, 5, 30]. This is indeed at the basis of the recently proposed 
ghost-free theory of massive gravity. Here however we are study- 
ing the theory with just a mass term. 
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With a standard analysis, discussed in App. A, we can 
recover from this Lagrangian the known result that one 
of these two degrees of freedom is a ghost while the other 
has the correct sign for the kinetic term. In particular, 
for (&i + b 2 )/b 1 b 2 < 0, the "healthy" mode is $ while 
the ghost is T. 6 Concerning the masses of the two scalar 
modes, as we show in App. A, the poles of the propagator 
are at 



m 



2 &i(6i+46 2 ) 2 



2(6i 



(77) 



In particular, setting b\ = l,fe 2 = — 1 — e the second pole, 
which corresponds to the ghost, is at 



3 + 4e 2 



(78) 



2. Approaching the FP point, b± + bi — > with b\ > 

Equation (78) shows that the ghost mass diverges as 
e — > 0. In this limit we can also easily see this directly 
from the Lagrangian, writing b\ = 1 and 6 2 = — 1 — e, 
and neglecting terms that vanish as e — > 0. Then eq. (76) 
becomes 



£ 2 , S calar($,r) = 12 



1 + 2e 
32e 



2t^2 



m 2 r$. 



(79) 



Note also that the term proportional to di&d 1 ® disap- 
peared, and therefore the ghost is simply T, see foot- 
note 6. Diagonalizing the kinetic term as in App. A, 



C 



2. scalar 



($,r) 



i 

48 1 



124$ 



12m 2 <I> 2 



1 

48 1 
-2e 



n 2 



12a $9*$ 



2t-.2 

m„L 

32e 9 



m 2 r$. 



(80) 



As e — > the mass term for T diverges, and T effectively 
disappears from the theory. This result is independent 
of whether we take the limit e — > + or e — > 0~ , i.e. of 
whether the mass term for T is tachyonic or not. This 
can be seen more clearly observing that the equation of 
motion derived from the Lagrangian density (79) taking 
the variation with respect to T is 



-e($ 



m 2 $) 



11) 



which in the limit e — > gives L = 0. Therefore at the FP 
point the ghost disappears and, taking the limit e — > in 
eq. (79), we remain with 



C 



2, scalar 



12 



$ 2 - & $5*$ - m 2 $ 2 

1 9 



(82) 



This Lagrangian density has the correct (i.e. non ghost- 
like) sign for the kinetic term, and the correct non- 
tachyonic sign for the mass term. The equation of motion 
for $ at the FP point is the massive Klein-Gordon (KG) 
equation, 



(□ -m 2 )$ = 0. 



(83) 



Thus the Bardecn "potential" $, that in the masslcss 
case is a physical but non-radiative degree of freedom, at 
the FP point becomes a radiative field which describes 
the helicity-0 component of the massive graviton. This 
is a result that was obtained already in 1966 in [24] by 
performing the elimination of the non-dynamical degrees 
of freedom in the ADM formalism, 7 and more recently in 
[26]. 

At the technical level it is remarkable that, at the FP 
point, the complicated Lagrangian density (76) collapses 
to the above simple KG form, particularly considering 
that $ is not a Lorentz scalar. We have indeed seen in 
eq. (55) that under boosts it mixes with S\ In partic- 
ular, the Lagrangian (82) is not Lorentz invariant. The 
Lorentz invariance of the theory only emerges when we 
recombine the scalar, vector and tensor sectors. However, 
the appearance of a KG equation is expected on physical 
grounds since, once we eliminate all the non-dynamical 
variables, the remaining physical degree of freedom must 
satisfy a Lorentz-invariant dispersion relation. 



3. Working directly at the FP point bi = 1, 62 = —1. 

It is instructive to derive the above results working 
directly at the FP point, rather than approaching it as 
a limit, since the structure of the equations is different. 
The variable 7 can still be eliminated using the variation 
with respect to 7, that now becomes 



4$ 

TO 2 



(84) 



The equation obtained from the variation with respect to 
ip now has a different structure: the term proportional 



For (61 + i>2 £>2 > the situation is more involved. Be- 
cause of the term c9;<I > c9 l <I> in eq. (76), as we discuss in App. A 
there is a critical value fc* = |k*| defined by (k t /m g ) 2 = 
(3/4) (b\ b'2 ) / (61 + 62 ) , such that are ghost-like the Fourier modes 
Tk with |k| < k„ and the Fourier modes ojfc<J>k+rk with k| > k„, 
where is denned in eq. (A3). However, even in this case, as 
bi + b'z — > we have fc* — > +00 and only the Fourier modes of T 
arc ghost-like. 



7 In [24] the scalar sector of the spatial metric was parametrized 



as hi 



{l/2)(8ij - 'V- 2 d i di)h T + 2didjh L and it was indeed 



found that the Hamiltonian has the KG form in terms of the 
variable h T and of its conjugate momentum. In our notation 
h T = -Acj> - (2/3)V 2 A = 4$. 
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to ip 2 in eq. (68) cancels, and ip enters in the Lagrangian 
£ 2 , SC aiar only linearly, 



2, scalar 



A.scalar + V> $ ~ 12m»$ - 2m 2 V 2 A] 



(85) 

where C 2 sca i ar is independent of ip. Thus ip is a Lagrange 
multiplier that imposes the condition 



3V 2 $ - 12m 2 $ - 2m 2 V 2 A = . 



(86) 



which is nothing but T = 0. The variation with respect 
to ip therefore eliminates the ghost. After eliminating 7 
through cq. (84) and using the condition T = we remain 
with a Lagrangian function of $ only which, as expected, 
gives back eq. (82). 



4- Interaction with matter sources at the FP point 

We now include the coupling to the energy-momentum 
tensor, to see how the elimination of the non-dynamical 
variables works in the presence of matter and what quan- 
tity acts as a source for the radiative field <£>. From 
eq. (31) the interaction Lagrangian in the scalar sector is 



Ant.scalar = «(3$P + * p) 



(87) 



where the factor of k comes from the fact that here the 
harmonic decomposition is performed on the metric per- 
turbation /i M „ defined by g^ v = r/^ + nh^, while in Sect. 
II we were using g^ v = rj^ + h^ v . Adding this interac- 
tion term to the Lagrangian in eq. (68), specializing to 
the FP point, and writing the Lagrangain in terms of the 
independent variables {7, ip, A, 4>}, we get 

(£ 2 +Ant) S calar = -12$ 2 + 45,$^$ - 8^^$ 

+8a t 7<9 l $ + 4$V 2 A 

+m 2 [12$ 2 + 4$V 2 A + <9. t 7d'7 - 12$V - 2^V 2 A] 
+ ~[6<f>P + 2ipp-2jp + SV 2 X}. (88) 

In the last line we used the conservation equation (26) 
and integration by parts to replace Xp — > Xp = AV 2 5 — > 
5V 2 A. The variations with respect to j,ip and A give, 
respectively, 

2to 2 7 ) = -up. (89) 
(90) 



8V 2 $ - 12m 2 g <P - 2m 2 V 2 A 



-K P . 
1 



2§ + 2m*<S>-m 2 g ip = ~2 kS - ( 91 ) 

Using again the conservation equation (26), p = V 2 5, 
we can rewrite eq. (89) as 



m 2 7 = 4$+-5. 



(92) 



We now eliminate 7 from the Lagrangian density 
using this equation, and V 2 A and ip using eqs. (90) and 



(91), respectively. We get £ scalar = £ 2 , sca iar + Ant.scaia 
with 



^scalar = -12(^$^$ + m 2 $ 2 ) 



and 



Cs 



int. scalar 



K$(2p + 3P - 35) 



(93) 



(94) 



Observe that all terms quadratic in the sources, which 
appear in the intermediate steps of the computation, in 
the end canceled. The equation of motion for $ in the 
presence of sources is therefore 



(□-m 2 )$ 



24 



{2p + 3P - 35) , ( [$] ~ mass) , (95) 



where we have recalled that in this section the variable 
$ was defined from the harmonic decomposition of the 
field h^y defined from g^ v = i]^ + nh^, so both h^ v and 
<3> have the canonical dimensions of mass. If instead we 
use the more common definition of $ from the harmonic 
decomposition of the field defined from g^ u = r/^ + 
(as we did in Sect. II), then both h^ v and $ are 
dimcnsionlcss, and eq. (95) becomes 



(□ - m 2 )* = 
v 9 ' 24 



-(2p + 3P-35) Q$]~l).(96) 
Using eq. (29), we can also rewrite this as 

([*]-!), (97) 



(□ - TO 2 )$ = — (p - V 2 cr) 



showing that, at the level of linearized theory, the 
hclicity-0 component of the massive graviton is sourced 
by the combination (p — V 2 cr). Observe that in con- 
trast, in the massless case, <I> was coupled only to p, and 
through a Poisson equation. 



5. The higher-derivative action for $ 

We now follow a different route and, for 61, 6 2 generic, 
we eliminate T and obtain a higher-derivative Lagrangian 
for <E>. From eq. (76) we see that the variation with re- 
spect to T gives 



16(6! + 6 2 ) 



61(61 



26 2 )m2 



m g bi$ 



61 



^V 2 $ 



(98) 

Thus, as we already remarked, even T can be elimi- 
nated algebraically, but the price to be paid is a higher- 
derivative action for We plug this expression for T 
into eq. (76) and, after some integration by parts, we find 
that all terms "miraculously" combine to form explicitly 
Lorentz-covariant structures, 



£ 2 ($) = ond^d^ + a 2 m 2 J>' 2 



03 



(99) 
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with 



Oil 



4(bi - 26 2 ) 
h + 2b 2 



Oi2 



4&i(bi + 46 2 
61 + 26 2 



a 3 



8(61+62) 
b 1 (b 1 +2b 2 ) 



(100) 



(101) 



At the FP point the higher-derivative term disappears 
and the result reduces to eq. (82), as it should. The 
term (l/mp(D$) 2 shows that, without the FP tuning, 
the theory becomes strongly coupled already at the scale 
m g . This is a result which is usually derived introducing 
a scalar Stiickclbcrg field, but which is quite transparent 
also in our formalism. 

Including also the interaction with the energy- 
momentum tensor, and repeating for 61,62 generic the 
computations done at the FP point in sect. Ill A 4, after 
long but straightforward algebra we get 



Ant(*) = K$ 



26 2 6l+ 462 A 

P + or — rr— o 



+ 



61 +26 2 

2(6l +b 2 )K 

61(61 + 26 2 )m 2 



61 + 26 2 

(p-S)n<z> 



+ 86 1 (6i + 26 2 )m 2 / P j ' 



(102) 



As it should, at the FP point the result reduces to 
eq. (94). Observe however that outside the FP point 
there is also a term quadratic in the sources, proportional 
to (p-S) 2 . 



B. Elimination of the non-dynamical fields. Vector 
and tensor sector 



A similar analysis can be performed in the vector sector 
of the theory, see also [25, 26]. With the addition of a 
generic mass term the quadratic part of the Lagrangian 
density in the vector sector reads 



£ 2 ,vcctor - diS^+mp! [ fop - jdtej&el ) , (103) 



where Si = (l/2)ej. We use /3j and ej as independent 
fields. We see that fa is a non-dynamical field, and the 
variation with respect to it gives 



2V 2 /? 1 - V 2 e l - 2m 2 b 1 f3 l = . 



(104) 



In contrast ej is dynamical, since Sj contains ej, and its 
variation gives 



V 2 e i -2/3 i + m 2 6ieM =0 



(105) 



Taking the time derivative of eq. (104) and substituting 
e' from eq. (106) we get the relation 



2 j 



2f3 l 



(107) 



which inserted back into eq. (106) gives 

(□ - m 2 G )e l = , (108) 

where mg = 61m 2 . Thus, the mass of the two vector 
modes described by the transverse vector e l is the same 
as the mass of the scalar mode. Clearly, e l describes the 
two components of the massive graviton with helicitics 
±1. 8 

For 61 < 0, however, these vector modes are not only 
tachyonic (since = 61m 2 becomes negative) but even 
ghost-like. The simplest way to understand this point 
is to express the Lagrangian in terms of the momentum 
modes, performing a spatial Fourier transform. From 
eq. (104), for the Fourier modes p\(i) and e^(t) we get 



4 



1 



2 k 2 + m 2 6i 



(109) 



This allows one to eliminate p\ from the action. Then, 
for the Lagrangian L = J d 3 xC one obtains [25] 



L 



2. vector 



m 2 6i 



d 3 k 



(2tt) 3 k 2 



61m 2 



(110) 



£i,k(k 2 + 6im 2 )4] 



If 61 > we have tuq 

-^2,vector(Ci) = " 



6lTO 2 > 0, 



d 3 k 



and 



(2tt) 3 



i,k c k 



-i,k 



k 2 + 
(k 2 



(111) 



; )4 



We see that the kinetic term has the good, non-ghostlike, 
sign and the mass term is non-tachyonic. For each given 
momentum mode e^, the equation of motion derived from 
the Lagrangian (111) is the same as that derived from 

(k 2 + m^e^ , which is just 



e i,k e k 



e i,kV 



a Lagrangian 

a massive Klein-Gordon Lagrangian, so we recover the 
massive KG equation (108). In contrast, if 61 < the 
sign of the kinetic term in the action is such that the 



modes of et with Ikl > 



-61 to 2 are ghost-like. 



This 



eliminates the possibility of a consistent theory at the 
sign- reversed of the FP point, 61 = — 1, 6 2 = 1. 

Observe also that, going back to coordinate space, the 
Lagrangian (111) becomes [26] 



L 



2, vector 



(112) 



With the boundary condition that the fields vanish at 
infinity, this is equivalent to 



e* - 2j8* + ro'&ie* = . 



(106) 



8 We adhere here to a common abuse of language. Of course, these 
modes become helicity eigenstates only in the massless limit. 
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The non-locality of this expression with respect to the 
spatial coordinates is ultimately a consequence of the fact 
that e l is a non-local function of the metric. However we 
have seen above that in the end, the equation of motion 
for e l is just a massive KG equation, so it is local. 

We next add the interaction term. Using eq. (31) we 
see that in the vector sector (with the present normaliza- 
tion of the fields) it is given by £j nt = — KpiS' 1 + (k /2)eiS l . 
Then eqs. (104) and (105) become 

2V 2 /3 l - V 2 e ; - 2m 2 G ? = -kS* , (113) 
V 2 (e* - 2$* + m 2 G e^ = kS 1 . (114) 

Using the conservation equation (30) and the usual 
boundary conditions that the fields vanish at infinity, 
eq. (114) becomes 

i* - 2p + m^e* = -<r\ (115) 

As in the sourceless case, taking the time derivative of 
eq. (113) and combining it with eq. (114) gives V 2 e l = 
2/3\ Plugging this into eq. (115) we get 

P = (H6) 

Finally, the tensor sector is simple, since we only have 
hJ 3 T . Using eq. (31), the Lagrangian density in the tensor 
sector is 

(£2 + £int) tensor = (H?) 

Again, b\ > gives a non-tachyonic mass term, while for 
61 < we have a tachyon. The equation of motion is 

(p-m&h*™ = (118) 

Comparing with eq. (77) we see that even when b± =/= I 
the scalar $, the vector e 1 and tensor modes h^' , that 
together make up the five components of the massive 
graviton, have the same mass, given by m G = 61m 2 . 
This is in agreement with the representation theory of 
the Poincare group, as it should for a theory linearized 
over flat space. 

IV. A HIDDEN GAUGE SYMMETRY 

The results of the previous section show that, out of the 
10 components of the metric, the four variables 7, ip and 
(3 l (where f3 z is transverse and therefore carries two de- 
grees of freedom) are non-physical and can be eliminated 
through their equations of motion. The five components 
of the massive graviton are described by {<&, e,, hj?}, 
with e 1 transverse and hj^ transverse-traceless. Finally, 
the ghost is given by a linear combination of V 2 A and 



and disappears at the FP point. In particular, when 
(61 + b2)/bib2 < 0, the ghost is given by the combination 
T defined in eq. (74). Alternatively, both scalar degrees 
of freedom can be described by the field $, through a 
higher-derivative action. 

We can describe this as a reduction process from 
the original action S[h^ v ], which is a functional of all 
the 10 components of the metric, to a reduced action 
SVcdl^, e»j hj?] which depends only on the physical de- 
grees of freedom (and includes a higher-derivative term 
proportional (D<I>) 2 if we are not at the FP point), 

S[V] (H9) 

This result, however, brings a surprise. Because of the 
mass term, the starting action Slh^] is not gauge invari- 
ant. Nevertheless, the variables $ and /i y ' TT are gauge 
invariant, while e 1 transforms as e — > ej— 2.Bj, see eq. (10). 
Therefore, none of these fields is sensitive to the transfor- 
mations parametrized by A and C , so the reduced action 
is trivially invariant under diffcomorphisms of the form 

£ =A, Ci = diC. (120) 

Observe that this is true independently of the form of 
the mass term. More precisely, if we perform a gauge 
transformation (120), the original action Slh^] trans- 
forms into a different quantity S"[/i M „; A, C]. However, 
as illustrated graphically in fig. 1, the reduction to the 
physical degrees of freedom performed with S'lh^A, C] 
(using of course its own non-dynamical equations of mo- 
tion) gives the same reduced action as that obtained from 

S'[h^;A, C] — > S rcd [$, e h hj?] , (121) 

independently of A and C. We therefore have a "hid- 
den" symmetry at the level of the original action: even 
if Slhftv] is not invariant under a gauge transformation 
of the form (120), still it is transformed to an action 
S'[/i M „; A, C] that describes the same physics. In this 
sense, eq. (120) is a symmetry of the original action 
S[h^]. 9 



To avoid misunderstanding, observe that we cannot use A and 
C to fix a gauge in the original theory, e.g. setting A = 7 = in 
the equations of motion, as we could do in the massless theory. 
In fact, using the Lagrangian £, the variables A and 7 are fixed 
in terms of <& by eqs. (84) and (86) (or by their generalizations 
outside the FP point), and cannot be set to zero. Using the 
Lagrangian £' they are instead given by eqs. (f25) and (126) 
below. In this case one might choose V 2 C and A such that 
A = 7 = 0. However, the resulting functions A and C would 
now depend on <E>, and would contribute to the action (125). 
Thus, the symmetry that we are discussing is different from a 
usual gauge symmetry, in which the action is invariant under 
the gauge transformation and the gauge freedom can be used to 
remove some degrees of freedom from the theory. In our case 
the original action is not invariant, but changes to a new action 
S'\hau', A, C]. The only sense in which these transformations are 
a symmetry of the theory is the one illustrated graphically in 
fig. 1. 
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sPv v ] 



A,C 



S'[h^ v ;A, C] 




S red [<£,8i,hV ] 




Figure 1: A graphic representation of the relations between 
5, S' and the reduced action. 



To illustrate this equivalence we consider the original 
action in the scalar sector, eq. (68), and we study how it 
changes under a transformation parametrized by A and 
C. We work with bi, b 2 generic and we use {$, i/), 7, A} as 
independent fields. We then transform these variables ac- 
cording to eqs. (8)-(10) and we find that the transformed 
Lagrangian density is 

4, 8calar [^V,7 J A;A,C]= (122) 

^2,scalar[$, lp, 7, M + A£ 2 ,scalar 0>, 1, A; A, C] , 

where £2. scalar is the original Lagrangian. Written in 
terms of the variables {<£, ip, 7, A, } it reads 

£a,Bcalar[*,V»7,A]= (123) 

-12$ 2 + 4di$<9 4 $ - 8^V9 4 $ + 8<%7<9*$ + 4$V 2 A 

- 2(6i + & 2 )^ 2 ~ 6(61 + 36 2 )$ 2 - 2(6i + 3& 2 )$V 2 A 

(61 + b 2 )(V 2 X) 2 + hdijd^ + 126 2 $V + 2& 2 0V 2 A 



(124) 



The extra term is 



.scalar 

[$,0,7,A;A,C] = 
m 2 i [4(6i + & 2 )V> - 126 2 $ - 26 2 V 2 A] - 2bim 2 g d l Ad t j 
+m 2 V 2 C [2(6i + 6 2 )V 2 A + 4(61 + 3fe 2 )$ - 2fe i7 - 46 2 0>] 
-2(6i + b 2 )m 2 g (A - V 2 C) 2 + hm 2 g d l (A - C)d t (A - C) . 

We can now take the variation of £2, scalar + A£ 2 sca i ar 
with respect to 7 and ip. Specializing for simplicity to 
the FP point we get, respectively, 



m-7 = 4$ + m 2 (A + C), 



m^V^A = 4V$-6m;$ + 2m;VC, 



(125) 
(126) 



Not surprisingly, these equations could have been ob- 
tained simply performing the replacements eqs. (8)-(9) 
directly into eqs. (84) and (86). Substituting these ex- 
pressions for 7 and V 2 A back into eqs. (123) and (125) 
(with b\ = 1,& 2 = —1) we find that, in £ 2jS caiar + 
A£ 2iSca i ar , all terms that depend on A and C (both linear 



and quadratic) cancel, and therefore the reduced action 
is the same as that found in the case A = C = 0. The 
same calculation can be performed for 61 + 6 2 7^ 0, and 
we have checked that again the extra terms cancel upon 
use of the equations of motion that eliminate 7, and 
V 2 A (or, equivalently, T). 

The symmetry that we have found takes a simpler form 
if, rather than using {A, C) as independent gauge func- 
tions, we introduce 



A = A-C. 



(127) 



and we use the pair (A, C) as independent variables. 
Then the residual gauge symmetry consists of linearized 
diffeomorphisms with 



H, = A5l + d,C. 



(128) 



In particular the transformation parametrized by C takes 
a nice covariant form, 



2d„d u C . 



(129) 



It is interesting to understand the origin of this hidden 
symmetry using this covariant form. Under h^ v — > h^ v — 
2d^d y C the term /i M!/ £ A " y,ptT h pa in eq. (65) is obviously 
invariant, and transforming the mass term we get 

4tV;<?] = \ [KvS^Kv - mKhh^hr + b 2 h 2 )] 

-2m 2 [(dvQd^hW + b 2 jf"h) + (61 + b 2 ){UC) 2 )} . 

(130) 

Observe that, at the FP point, the term quadratic in 
C drops. This corresponds to the well known fact that, 
when hu V takes the form hau — d^d u 4> for some function 
</>, the FP mass term is a total derivative. The generaliza- 
tion of this property to terms of cubic and higher order in 
d^dycj) gives rise to the Galileon family of operators [31]. 

In this covariant formulation a simple but not to- 
tally correct way of showing that, upon use of the non- 
dynamical equations of motion, the extra term in the 
Lagrangian vanishes, is as follows. Define 



f ^d^bw + bzrrh). 



(131) 



We see from eq. (67) that, if we use the equations of 
motion of the starting Lagrangian C 2l we have f = 0. 
However, as we show in app. B, to derive this condition 
it is not sufficient to use the non-dynamical equations 
of motion. Rather, we also need the equation of motion 
for e l , which is a dynamical variable. Still, using only 
non-dynamical equations of motion, we can show that 



and therefore also 



0. 



0. 



duj" = 0. 



(132) 



(133) 



Thus, as a consequence of the equations that eliminate 
the non-dynamical degrees of freedom, the term (d t ,C)j l/ 
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in eq. (130) vanishes upon integration by parts. The 
term (DC) 2 is anyhow irrelevant, even away from the FP 
point, since it is decoupled from the other fields, and we 
can seemingly conclude to the invariance of the reduced 
theory. 

The reason why this argument is not totally correct 
is that, to perform the reduction of £ to the physical 
degrees of freedom, we must use the equations of mo- 
tion of £' , rather than the equations of motion of C. 
Thus, the correct derivation is the one that we have given 
above, in which we have have eliminated from £' the non- 
dynamical variables, using their equations of motion de- 
rived from Lagrangian C itself. In the end, however, the 
only difference between using the equations of motion of 
£' or that of C is that with the former procedure (which 
in principle is the correct one) both the term linear in 
(A, C) and the terms quadratic in (A, C) cancel exactly, 
as we checked explicitly. In contrast, with the (a pri- 
ori incorrect) procedure of using the equations of motion 
obtained varying C to eliminate the non-dynamical vari- 
ables from the terms linear in A, C still cancel while 
those quadratic in A, C do not. Since these quadratic 
term are decoupled from the physical fields {<!>, e^, /i^ T } 
they are anyhow irrelevant, so in the end even the simpler 
procedure gives the correct answer. 



The overall Stiickclbcrg replacement is therefore 

1 2 
huv^hu V -\ {9a7t v + d v TCu) -\ ^8^71. (140) 

m g nig 

Writing £ M as in eq. (7), the transformation of the Stiick- 
clbcrg fields under linearized diffeomorphisms, eq. (139), 
becomes 

7t n + m g (A - C) , Tti ->• 7tj + m g Bi , (141) 



7W7T+TO 2 C, (142) 

while in terms of the harmonic variables the Stiickclbcrg 
transformation (140) reads 



1p -> 1p + — 7T0 H yfC. 



-> 0- 



n>7, 



V 2 7t, (143) 



3m 2 



7^7 



-7T 



:Tt. 



A — s- A + 



(144) 



(145) 



V. HIDDEN SYMMETRY IN THE 
STUCKELBERG FORMALISM 

An equivalent and instructive derivation of the exis- 
tence of this hidden symmetry can be given by making 
use of the Stiickclbcrg formalism. We first introduce 
as usual a Stiickelberg vector field n^, performing the 
replacement 

hpv -> h^v -\ {d^n v + d v n,j) (134) 

m g 

in the action for linearized massive gravity. The lin- 
earized massive theory is then trivially invariant under 
the combined transformation 



, + . (136) 



7t„ 



We next introduce the Stiickelberg scalar n, replacing 
further 



1 A 

Ttjn ->• Tlfi H d^Tt 



(137) 



The theory then also acquires a 17(1) gauge symmetry, 

Ttfi Y TC^ d^6 , 7T-> n + m g e, (138) 
while eq. (136) becomes 



nig nig 



(139) 



We can now see how the hidden symmetry emerges in 
this formalism. When we perform the replacement (143)- 
(145) into the Lagrangian eq. (68) we obtain a new La- 
grangian that depends both on the metric and on the 
Stiickelberg fields. We then write down the equations of 
motion of this Lagrangian with respect to 7, ifi and V 2 A 
and use them to eliminate these variables from the the- 
ory. Obviously, since the replacements (143) and (144) 
are patterned after the gauge transformation of the scalar 
fields, this is formally the same computation already dis- 
cussed in sect. IV, with n playing the role of — m 2 C and 
7to of — m g (A— C). Thus the result is that, after elimina- 
tion of the variables ip, 7 and V 2 A, the Stiickelberg fields 
7t and Tto disappear from the action, just as all terms in- 
volving A and C canceled in the computation discussed 
in sect. IV. In conclusion, after reduction to the five 
physical degrees of freedom of the massive theory, the 
Stiickelberg fields tcq and 7T drop from the action, and 
are no longer needed for restoring the scalar part of the 
gauge symmetry. In contrast, 7tj remains in the vector 
sector, where it is needed to obtain invariance under the 
transformation parametrized by 

Observe also that, if we use the pair (^4, C) as inde- 
pendent variables, the Stiickelberg fields Tto and n take 
care separately of the transformations parametrized by 
A and C, 



7t 



Tto + m g A , 



7t 



n + nigC . 



(146) 



In this formulation, where the transformation 



parametrized by C is simply hn V 



we immediately recognize that the hidden symmetry 
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parametrized by C is nothing but the U(l) gauge sym- 
metry (138). Indeed, from eq. (134), the transformation 
7t^ — y 71^ — d^B induces on the metric the transformation 
h^v — > h^ v — 2d fl d„C with C = 9/m g . The reason 
for performing a further Stuckelberg transformation, 
introducing the Stiickclberg scalar n, was to cancel this 
variation. However we have seen that, after reduction 
to the physical degrees of freedom, the transformation 
h^u — > h/ju/ — 2d fi d„C becomes a symmetry of the theory. 
Thus, after performing the Stuckelberg transformation 
(134) and eliminating the non-physical components of 
the metric, the U(l) gauge symmetry is already there, 
without the need of introducing further the Stuckelberg 
scalar TC. 

VI. (3+1) DECOMPOSITION AND BARDEEN'S 
VARIABLE IN ELECTRODYNAMICS 

It is interesting to compare the above results with a 
similar analysis performed both in massless and in mas- 
sive electrodynamics. We begin with the massless theory 
and apply the harmonic decomposition to the electro- 
magnetic field, writing the gauge field as 

A = ip , A i =v i + d i X, (147) 

where Vi is a transverse vector, diV 1 = 0. Under a gauge 
transformation — > A^ — d^O we have 

ip^-ip-0, A ^ A - (9 , (148) 

while Vi is invariant (similarly to the fact that in lin- 
earized gravity hj^ is invariant). We can then define a 
"Bardcen variable" 

$> = ip-X, (149) 

which is gauge invariant. The four degrees of freedom 
of A^ have therefore been rearranged into three gauge- 
invariant degrees of freedom {^, Uj}, and one pure-gauge 
degree of freedom, that can be taken to be A or -0. In- 
verting eq. (147) we get 

\ = V- 2 (d l A l ), v i = A i -d i V- 2 (d j A j ). (150) 

Thus, as in linearized gravity, these harmonic variables 
are non-local function of the gauge field A^ and the non- 
locality disappears in the gauge diA 1 = 0. We now couple 
the electromagnetic field to an external current j^, and 
by analogy with eqs. (22) and (23) we write 

jo = P , ji = Si + d l S , (151) 

with diS 1 = 0. Current conservation implies V 2 5 = p. 
After some integrations by parts and use of this conserva- 
tion equation, the Lagrangian density for the electromag- 
netic field coupled to an external current can be written 

as 

= i [flsttfl** - fy^dV] + (js9 - Sirf) . (152) 



As expected, the Lagrangian density of the massless the- 
ory depends only on the gauge-invariant variables $ and 
Vi. The variations with respect to $ and Vi give, respec- 
tively, 

V 2 f = p, □«* = £*. (153) 

This shows that v % describes the two radiative degrees of 
freedom of the massless photon, while ^ is a physical but 
non-radiative degree of freedom. 

We next consider the massive theory defined by the 
Proca Lagrangian density 

C = --F^ v - im^A" . (154) 

Expressing it in harmonic variables we get C = £ sca iar + 

^vector with 

Scalar = \ [$#0** + m* - d t Xd* A)] , (155) 

^vector = --[d lt v i & l v i +m*ViV i ] . (156) 

In the vector sector the variation with respect to Vi gives 
the expected massive KG equation, 

(□-m^y = 0. (157) 

The Lagrangian in the scalar sector now depends sep- 
arately on ip and A, rather than only on the gauge- 
invariant combination \&. A difference with respect to 
linearized gravity with the FP mass term is that now the 
variable ip enters quadratically. We use A and ip as in- 
dependent variables. Clearly ip is non dynamical, since 
(recalling that $ = ip — A) it enters the action without 
time derivatives. The variations with respect to ip and A 
give, respectively, 

V 2 ip - m 2 ip = V 2 A , (158) 

A + ?n 2 A = ip . (159) 

Taking the time derivative of eq. (158) and replacing ip 
from eq. (159) we get 

(□-m 2 )A = 0. (160) 

This shows that the radiative degree of freedom in the 
scalar sector is A, rather than the gauge- invariant vari- 
able \& (an obvious result, since we see from eq. (147) that 
A corresponds to the longitudinal polarization). The vari- 
able W is instead determined by A through the equation 
(V 2 — m 2 )^ = miX, which is obtained writing ip = \S r +A 
in eq. (158). Since A transforms as A — > A — 0, even af- 
ter eliminating the non-dynamical variable ip no gauge 
symmetry appears in the Proca formulation of massive 
electrodynamics. The "hidden" gauge symmetry is really 
a peculiarity of linearized massive gravity. 
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It is also instructive to see how the Lorentz invariance 
of massive electrodynamics is recovered using the har- 
monic variables. Under an infinitesimal boost we have 
A Q -> A Q + uj Qi A l and A 1 -> A 1 + uo l A a . Then, proceed- 
ing as we did in sect. II B 2, we find that under boosts 

Syj = uoHvt + diX), (161) 
SX = Lo Q l V- 2 [d t y + v,], (162) 
5vi = -Lu ^V- 2 [{didj - %V 2 )* + diVj] ,(163) 

and therefore (recalling that 5X = 6(doX) = uirfdiX + 
d (SX)) 



5^ = ujo l V (Ovi-diV 



(164) 



Observe that, just as in the gravitational case, the gauge- 
invariant variables and u,: transform among themselves 
under Lorentz transformations, but their transformation 
involve the non-local operator V~ 2 . 

Using 6(diV) = 6{d t )^ + d t S^ = oj 0i ^ + d t S^ we find 



(165) 



(plus a term dido [(l/2)wo 4l 3/ 2 ] which is a total derivative). 
Similarly, neglecting total derivatives, we get 



5 (d^d^v 1 ) = -2w i *nw. t . 



(166) 



and 



8{ip 2 - d t Xd l X) = S(v t v l ) = 2uj Q i ^v i . (167) 
Therefore, modulo total derivatives, 

^scalar = -^vector = -UJ 0i ^ (O - m 2 g )v % . (168) 

Thus the variation of £ S caiar cancels the variation of 
>Cvector and the total action is Lorcntz-invariant, as it 
should. Observe that Scalar and ^vector are not sepa- 
rately Lorentz invariant (unless we impose the equation 
of motion for v l ). 



VII. RADIATIVE DEGREES OF FREEDOM 
AND HIDDEN GAUGE SYMMETRY IN 
DE SITTER SPACE 

In this section we consider massive gravity linearized 
over a de Sitter background. We start from the Einstein- 
Hilbcrt action with a cosmological constant supple- 
mented with a FP mass term, and we work for generality 
in d spatial dimensions. Then 



S = / d d+1 x- 



(ij_2A)--£ {h^h^-h 2 
K 2 v 



(169) 

where is defined from g M „ = g^ v + nh^ LU . We write 
the background metric in the form g^ v = (— l,a (t)<5,j), 



with H = a I a = const. The linearization of the action 
gives [11] 



S 2 = I d d+1 



x v — g 



+ dfl { h ^ V \ k2 ) it ^ 

(170) 

where D M is the covariant derivative with respect to the 
background metric g^ u , and R is the Ricci scalar of the 
background. In de Sitter space with d spatial dimen- 
sions R = d(d+ 1)H 2 , and the Einstein equations for the 
background fix it in terms of the cosmological constant, 
R = 2[(d + l)/(d - 1)]A, so H 2 = 2A/[d(d - 1)]. 

As is well known, if m g and H satisfy the Higuchi 
bound [32] m 2 > (d — 1)H 2 this theory has no ghost 
and five physical degrees of freedom: one in the scalar 
sector, two in the vector sector and two in the tensor 
sector. For m 2 < (d — 1)H 2 the degree of freedom in the 
scalar sector becomes a ghost, while at the special point 
m 2 = (d—l)H 2 the scalar disappears and we remain with 
just four propagating degrees of freedom. The theory at 
the special point m 2 = (d — l)H 2 is known as "partially 
massless" , and acquires an extra gauge symmetry of the 
form [33-35] 



5h 



DnDvCt 



1 



(171) 



parametrized by a function a(x). It is interesting to re- 
cover these results using the (3 + 1) decomposition, and 
see if the hidden symmetry found in the linearization over 
Minkowski persists. We limit ourselves to the scalar sec- 
tor of the theory. We parametrize the metric fluctuations 
in the scalar sector as 



/i o = 2ip , h i = di-f . 



hi. 



-2a 2 



did. 



X. 



(172) 



(173) 



We will use the sum over repeated spatial lower indices 
to denote contractions performed with the flat Minkowski 
metric, while contractions between upper and lower in- 
dices are done with the background FRW metric, so 



dtd t 



S^didj, d l d l = a- 2 5»didj, and V 2 



didi 



still denotes the flat-space Laplacian. The factor 1/d in 
cq. (173) ensures that the operator acting on A is trace- 
less in d spatial dimensions. 10 The massless theory is 
invariant under linearized diffcomorphisms 



huu -> h 



(174) 



10 The functions E and B commonly used in the literature arc re- 
lated to 7 and A by 7 = aB and A = 2a? E. 
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In the scalar sector we write again £0 = A and & = diC. 
Then the gauge transformations of the scalar functions 
are 

, 1 



ip -> ip - A , (j>-> 
X^\-2C, j^j-A-C + 2HC. (176) 



-APC-HA, (175) 
a 



In three-dimensional Minkowski space these transforma- 
tions reduce to eqs. (8) and (9), as they should. The 
Bardcen variables in FRW read 

$ = -4> - ^d t d l \ + Hj - ^(\ - 2H\) , (177) 

* = tp-j+^(X-2HX). (178) 

We find useful to rewrite them as 

$ = $o + i?C, * = V-C, (179) 

where 

^ = -^-^d t d l X (180) 

is the Bardeen variable in Minkowski space and 

C = l-\\ + H\. (181) 

Observe that under gauge transformations £ — > £ — A. 
From eqs. (175) and (176) it is straightforward to show 
that $ and \& are invariant under linearized diffeomor- 
phisms. Using eqs. (172)-(173) we find that the action 
in the scalar sector is 11 



S 2 , 8 calar = 2 ( d ~ 1) / ^ [-<*(* + # *) 

+{d - 2)di<S>d l <S> - 28^8^ 



+ - 



d d+1 xa d 



4d(d - l)(f 2 + 8d(/)i(j 



+2d i ^d i -i-^—^{d i d i Xf \ . (182) 



Setting a(t) = 1 and d = 3, we recover the flat-space re- 
sult (68). We now eliminate the non-dynamical degrees 



of freedom from the action and identify the variable de- 
scribing the radiative degree of freedom. 12 As indepen- 
dent variables we use the set {$0, tp, 7, A}, which is the 
same that we used in the Minkowski space. Note that in 
de Sitter $ contains A and therefore the change of vari- 
ables from <j> to 3> would not be legitimate, see footnote 4. 

We now write the action (182) using these variables. It 
is useful to rewrite the second and third term using the 
relation 



J dta d [{d - 2)c\$<9 4 $ - 2^$^*] 



H V 



(183) 



which is proved writing 2a di<&d l <& = 2a 2 di$d l $ 



a d 2 do(di§di$) and integrating by parts. We also ob- 
serve that 4> + = <3>o + Hip. Then we get 

S 2 ,scaiax = 2(d-1) J d d+1 xa d 
2 2 



x 



d ($o + fty) - -d % ($ + Hi) d l ($ + m 
+2(d-l) / d d+1 xa d {d i d i X) (184) 

d-1 



($0 + + dH ($ + Hip) + m 2 $ - 



m 2 g / d d+1 xa d [2d(d - - 4dip& + ^7^7] 



where we have collected the terms proportional to did l X 
since it is a Lagrange multiplier. The equation of motion 
with respect to 7 is 



2(d- 1) 



$0 + Hijj 



(185) 



11 We have checked that for m g = eq. (182) agrees with the 
action given in eq. (61) of Rubakov and Tinyakov (RT) [25] , after 
performing the change of notation ip = fwTt 7 = a B, A = 2a? E 
and </> + (l/2d)d{d' l \ = V>RT and going from cosmic to conformal 
time. Observe however that our expression given in terms of 
Bardecn's variables is much simpler. The mass term examined 
in [25] is instead different from the one that we are considering, 
since even for a de Sitter background it is constructed writing 
the FP mass term for the variable g M „ — r^i,, rather than for 
9t±v — H/ns, and in this case the theory has a ghost, as shown in 
[25]. 



So, just as in the flat-space case, 7 is a non-dynamical 
variable that can be eliminated using its own equation of 



12 In App. C wc compare our discussion with the results presented 
in [27], where the equations of motion of the theory have been 
studied using Bardeen's variables, and a Klein-Gordon equation 
is obtained for the combination a _2 (<I> + 'J'). As we discuss in 
App. C, this however does not mean that a~ 2 ($ + ty) is the 
radiative degree of freedom of the theory in the scalar sector. 
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motion, and wc obtain 



^scalar = 2(d - 1) J d d+1 X CL d 



-d ( $ + ff^J - ZH^d^od 1 ($o + Hip 



+2(d-l) / d d+l xa d (d l d i X) ($ + 



+d# (*o + H^j + m^ - 
J d d+1 x a d [2d(d - l)$g - 4^$o] ■ 



We next rescale our variables 

{i>,$ ,\} ^ a-( d -V{i>,$ ,\} , 
define an effective mass 

M 2 =m 2 g - (d- l)H 2 , 

and use 

X = i>-(d- 1)$ 
as an independent variable instead of ip. We get 



(186) 



(187) 



(188) 



(189) 



.5' 



2,scala 



2(d-l) J d d+l xa- d+2 



-d[<S>o + H X ) -2H- 1 d l <P d l (<i> + Hx 
_2(d-l)^ + ff \ ^ / $Q + ^ 



M 2 

d-r 

2d 



+2(d-l) / d d+1 a;a- d+2 (5 i a i A) 



+2(d- 1) / d d+i a;a 



-dml<f> 2 



(190) 



For M ^ we now trade in the action 3> for the new 
variable 51 defined by 



n(t,x) = $ (t,x) + jrQT dt' x (^,x) 

(d- 



Af 2 



$o (to , x) + H $ (to , x) + #x(to , x) 



(191) 



where t is the time at which the initial conditions are 
given. The first two terms are chosen so that 



n = $ + H X . 



(192) 



The last term in brackets is time-independent and is 
a functional of the initial data $ (t ,x), 4> (ioj x ) and 
x(fo,x). The reason for this specific choice of the time- 
independent term will become clear below. Note that 
the initial data {$p(*o), $o(*o)} and {x(*o), x(*o)} fully 
determine {£l(to): ^(^o)}; 



mi (d - 1)H , 



fi(<o) = $o(io)+#x(io) 



(193) 



(where, for notational simplicity, here and in the follow- 
ing we do not write explicitly the x dependence) . Invert- 
ing we get 



M 2 (d-l)H ■ N 

*o(*o) = — n *o -- — r-n(*o), 
#o(fo) = n(t ) - H X (t ) . 



(194) 



The occurrences of $o other than in the combination 
(192) bring in terms that are non-local in time. Using 
eq. (194) we get in fact 



$o(t) = fi(t) - H di'x(f')) 



f2(to) + fm(t )) • (195) 



(d-l).ff 



However, in these new variables the constraint imposed 
by the Lagrange multiplier did l X reads 



d- 1 



(196) 



so integrating out A and eliminating x from the action 
makes eq. (195) local in time, 



$ = 



(d-l)ff 



o. 



(197) 



The time- independent term in the definition of O, 
eq. (191), was precisely chosen so that time-independent 
part in this expression vanishes. The action now depends 
only on f2 and, after some integrations by parts and re- 
arrangements we find 



2, scalar 



-2d(d- 1) 



(198) 



M 2 



d d+1 xa- d+2 [d^Qd^n + M 2 n 2 ] 



We now rescale back 



n 



-d-l 



n. 



(199) 



so that in the action we reconstruct the volume form 
y/=9 = a d , 

m 2 r 

2,scalar = -2d(d~l) — / d** 1 X 0* [d^fl + m^Q 21 



S- 



(200) 
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This also has the effect of replacing the effective mass 
M 2 by m 2 inside the square brackets. In terms of the 
original variables, the dynamical mode is 

fi(t) - *o(t) 

+a- ( - d - 1 '>{t)H f At'a d - 1 {t') [%b{t') - (d- l)$o(*')] 
(d-l)if / a(t) x 



M 2 V«(*o) 

x ($ (io) + #$o(*o) + Wo)) ■ (201) 
After an integration by parts, we can rewrite it as 

a(t) 

( > » (202) 



a(t 

tt($ o((o)+iJ#o) ) + | $o(to) 



Equations (200) and (203) are the main result of this 
section, and nicely display in a compact manner a num- 
ber of known features of massive gravity in de Sitter 
space. First of all we see that for M 2 > 0, i.e. for 
to 2 > (d — l)H 2 , the kinetic term of the scalar field 
has the "good" non-ghostlike sign and the mass term 
is non-tachyonic. For M 2 < we have a ghost instead 
(and the mass term in eq. (200) becomes tachyonic). For 
M = the action vanishes and the radiative degree of 
freedom in the scalar sector disappears. 13 Observe also 
that we smoothly retrieve the flat space-time result in 
the H — > limit. The action (200) also agrees with that 
found in [36] using a different route, namely embedding 
.D-dimcnsional dc Sitter space in (D + l)-dimensional 
Minkowski space and using the Stiickclbcrg formalism 
to isolate the helicity-0 mode. 

To make contact with the well-known result given by 
Deser and Waldron [35], for M 2 > we can introduce a 
new variable 



q a (t) = 2 



M 2 



1/2 



a d / 2 n(t) (203) 



and the action (200) can be rewritten as 



S 2 



1 



scalar 



d d +h 



d 2 H 2 

d^qod^qo + ( to 2 — j g5 



(204) 



13 Actually, the manipulations leading to eq. (200) were performed 
assuming M ^ 0. However, for M = 0, one can already eliminate 
X using the constraint imposed by the Lagrange multiplier A 
in eq. (190). This gives H\ = — ($o + H$o). Replacing this 
expression for \ into the action gives 52, scalar = 0, so in eq. (200) 
the limit M 2 — > is smooth. 



which, specialized to d — 3, reproduces the result of 
[35]. Note that using the variable qo the volume form 
y/— g = a d has been eliminated in favor of an additional 
contribution to the mass term. 

Finally, having found the expression for the field that 
describes the radiative degree of freedom in the scalar sec- 
tor, we can ask whether the gauge symmetry in the scalar 
sector is still preserved, as in Minkowski space. Equa- 
tion (203) shows explicitly that the dynamical variable 
Q is invariant under gauge transformations parametrized 
by C since, under this transformation, $o(i) and ip(t) 
are invariant. In contrast, under the transformation 
parametrized by A we have 



$o -» $o + HA , 



^> ip — A. 



(205) 



Therefore, the terms in f2(£) containing the combination 
$o(*) + Hijj(t) are invariant, while the term $o(*o) has a 
variation determined by A(to), so 



0(i,x) ft(t,x) 



a(t ) 



< d -V m 2 H 



(206) 



We see that the function f2(i, x) has a variation that de- 
pends only on the value of the function A(t) at the ini- 
tial time to, rather than on the whole function A(t). We 
therefore still get an invariance if we restrict to gauge 
functions A(t) that vanish at the initial time chosen 
for assigning the initial conditions. Until now we have 
worked with to arbitrary. A natural choice, however, is 
to chose to = — oo, and to restrict to gauge functions 
A(t) that vanish as t — > — oo. With this restriction on 
the initial value of the function A(t), the variable is 
invariant under the gauge transformations parametrized 
by both C and by A, and therefore the "hidden sym- 
metry" that we found in massive gravity linearized over 
Minkowski space persists in massive gravity linearized 
over dc Sitter space. 



VIII. CONCLUSIONS AND SUMMARY 

We conclude by summarizing the main results and 
equations of this rather long paper. The use of the 
(3 + 1) decomposition of the metric and of Bardeen's 
variables provides a valuable tool for understanding var- 
ious aspects of massive gravity. This formalism allows to 
identify and eliminate the non-dynamical degrees of free- 
dom from the action, working directly at the Lagrangian 
rather than Hamiltonian level, and to write down a re- 
duced action for the radiative degrees of freedom. We 
have shown how to carry out this elimination procedure 
both in Minkowski and in de Sitter space. For massive 
gravity linearized over Minkowski space (in d spatial di- 
mensions) we found that, in the scalar sector, the vari- 
able that describes the radiative degree of freedom of the 
massive graviton is the flat-space Bardeen variable 



1 

2d 



V 2 X, 



(207) 



19 



with cj), A defined by the harmonic decomposition of the 
metric, eq. (3). (In sects. II- IV we denoted this variable 
simply as $. Here we reserve the notation $ for the 
curved-space Bardeen potential). In terms of $07 f° r a 
mass term that has the Fierz-Pauli form and in generic 
d spatial dimensions, the reduced theory in the scalar 
sector is simply described by a Klein-Gordon action, 

^scalar = ~^{d ~ 1) / d d+1 X (5^*0^*0 + ™g$o) . 

(208) 

and therefore <f>o satisfies a KG equation, 



(□-m?)$o = 0. 



(209) 



a result that was already found using the ADM formalism 
in [24] (see also [26]). In the presence of external matter, 
we found that this KG equation is sourced by the combi- 
nation (p — V 2 ct), see eq. (97). This result is interesting 
and somewhat unexpected for various reasons. First of 
all, in the massless case the Bardeen variables describe a 
physical but non-radiative degree of freedom, as we re- 
called in sect. II B. It is surprising to see that, when we 
switch on a mass term, it is precisely the combination $o 
that describes the radiative degree of freedom. The real 
surprise, however, comes from the fact that a variable 
such as $o is gauge-invariant under linearized diffcomor- 
phisms. Thus, despite the fact that the mass term breaks 
the gauge invariance of the theory, after elimination of 
the non-dynamical variables the scalar sector of massive 
gravity is still gauge invariant, in the sense discussed in 
sect. IV. Writing the gauge transformation in the form 
h^v -> h^ v - (d^ v + d v £n), with £ = A, & = B l + d t C, 
the reduced theory is invariant under the gauge transfor- 
mations parametrized by the scalar functions A and C. 
The crucial point for the existence of this symmetry is the 
fact that, after eliminating the non-dynamical degrees 
of freedom, the scalar sector can be written uniquely in 
terms of <J> (possibly including higher-derivative terms, 
if we are not at the FP point). This was not obvious 
a priori: in principle, one could have remained with a 
different non-gauge invariant field. This is indeed what 
happens in the vector sector, where the propagating field 
is e l , rather than the gauge-invariant combination E*. 
Thus, in the vector sector, the symmetry parametrized 
by Bi is broken. Similarly, in massive electrodynamics in 
the scalar sector survives the longitudinal mode, which 
is not gauge- invariant. Thus, the existence of a hidden 
symmetry is a peculiar and non-trivial property of the 
scalar sector of linearized massive gravity. 

It is also interesting to explore the structure of the 
theory in these variables when the mass term deviates 
from the FP form. In this case the scalar sector can be 
described by two fields ($0j T) and is governed by the La- 
grangian density (76), which is second-order in the time 
derivatives. When the coefficients b\, 62 that parametrize 
the mass term satisfy (b\ +62)76162 < 0, the combination 
r given in terms of the metric in eq. (74) is the ghost, 
while $0 is the healthy mode. Otherwise the ghost is ob- 
tained from a mixture of the Fourier modes of $0 and T. 



Alternatively we can integrate out even T, at the price of 
a higher-derivative action. It is remarkable to see that in 
this case the complicated Lagrangian (76) collapses to a 
simple covariant form for $0 



£ 2 ,scalar = Ou^O^C 



-a 2 m 2 g $l 







(210) 



with coefficients ai, 012, «3 given (for d = 3) in eqs. (100) 
and (101), and in particular a.3 = at the FP point. The 
corresponding source term is given in eq. (102). The ap- 
pearance of such explicitly Lorentz-covariant structures 
is required by the fact that, in the end, the variables 
that describe the helicity-0 modes must have a Lorentz- 
invariant dispersion relation. Observe however that $0 
is not a scalar field under Lorentz transformations, but is 
scalar only under spatial rotations, and the Lorentz in- 
variance of the theory only emerges combining the scalar, 
vector and tensor sectors. 

Finally, we have explored massive gravity linearized 
over de Sitter using this formalism. In this case the elim- 
ination of the non-dynamical variables and the identifica- 
tion of the radiative degree of freedom in the scalar sector 
proved to be more subtle. The radiative degree of free- 
dom is not given by the curved-space Bardeen variable $ 
defined in eq. (177) (nor by the combination + de- 
spite the fact that it satisfies a KG equation, see app. C). 
Instead, it is described by the field £1 given in eq. (203), 
which has the form 

a d - 1 (t)n(t,x)= f dt'a d - 1 (^)($o+^)(i',x)+/(i ,x), 

Jto 

(211) 

where <£>o is still the flat-space Bardeen variable, 2ip = 
hoo, and /(io,x) is a time- independent function which is 
fixed by the initial conditions imposed at the initial time 
to- This expression is non-local in time. However, in the 
limit H — > 0, ft becomes local and smoothly reduces to 
$0 • In term of this variable, and for a FP mass term, 
the reduced theory in the scalar sector is again described 
by a simple KG action, 



S- 



2. scalar 



-2d(d- 1)- 



-g [d^ftd^n + m 2 g Q 21 



(212) 



This expression nicely summarizes a number of known 
facts about massive gravity in de Sitter: (1) if m 2 g > 
(d — 1)H 2 (the Higuchi bound) we have a single healthy 
degree of freedom in the scalar sector with the correct 
sign for the kinetic term and a non-tachyonic mass term 
(we always assume m 2 > 0). (2) for m 2 < (d — l)-ff 2 , 
the scalar degree of freedom becomes a ghost, and (3) 
at the special point m 2 = (d — l)H 2 , corresponding to 
the so-called partially massless case, the scalar degree of 
freedom disappears from the spectrum of the theory. 

Furthermore, the variable f2 is invariant under scalar 
gauge transformations parametrized by A and C (with 



20 



the extra condition that we restrict to gauge functions 
A{t) that vanish at the initial time to, which can be cho- 
sen to be to = — oo), and therefore the "hidden symme- 
try" that we have found in Minkowski space persists in 
de Sitter. 

Our analysis has always been performed without in- 
cluding the non-linearities of the theory. It would be 
interesting to generalize it to the non-linear ghost-free 
de Rham-Gabadadzc-Tolley model of massive gravity 
[4, 5]. Work on this is currently in progress. 
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Appendix A: Diagonalization of the kinetic term in 
the scalar sector 



In this appendix we study the diagonalization of the 
kinetic term of the action (76). To study the existence 
of ghosts we consider at first only the terms that depend 
on time derivatives, i.e. we consider the Lagrangian 



L 



d A x 



i2$ 2 + r$ 



16(6i 



6i& 2 m 2 



(Al) 



(so that the corresponding action is S = J dtL). It 
is convenient to perform a spatial Fourier transform of 
<&(t,x) and r(t, x) and work with the modes $k(i) and 

r k (t), so 



L = 



(2^)3 



where 



a k = 24 - 



32(6i 



(A2) 



(A3) 



Since <&(i, x) and r(i,x) are real, the momentum modes 
satisfy = <&_k and r£ = T_ k - The corresponding 
conjugate momenta are II$ k = SL/S&k = a-k^ + 
and IIr k = SL/ST^ = $£, and the Hamiltonian is 



H = 



d 3 k 

d 3 k 

(2tt) 3 



n$ k $k + nr k r k j -l 



(A4) 



so H = L. This Hamiltonian can be diagonalized writing 
it as 



H 



d 3 k 



2a k 



a k <S>i 



1 

2afe 



(A5) 

This shows that, independently of the sign of ak, the 
Lagrangian (A2) always has a ghost [25]. In particular, 



if ak > the ghost is Tk while for ak < are ghost-like 
the Fourier modes Tk + afc^k- We see from eq. (A3) that, 
if 



61 + 62 
61&2 



<0, 



(A6) 



ak is positive for all k, so the ghost is T(t,x). For (61 + 
b 2 )jb\b 2 > the situation is instead quite peculiar: there 
is a critical value fc* = |k*| defined by ak, = 0, i.e. 



m, 



3 h b 2 



4 61 



(A7) 



such that are ghost-like the modes Tk with |k| < k„ and 
the modes «fe$k + Tk with |k| > fc*. Thus, in the two- 
parameter space (61, 62) that parametrizes the mass term, 
there is a region where are ghost-like the Fourier modes 
of a combination of metric components with momentum 
k smaller than a critical value fc* , and the Fourier modes 
of a different combination of metric components with k > 
k*, a situation which is quite unusual. 

For comparison with the masses of the vector and ten- 
sor modes, it is interesting to compute the masses of the 
two scalar modes for b\ and 62 generic. This can be conve- 
niently done specializing to field configurations indepen- 
dent of x, i.e. to the Fourier modes with k = 0. We also 
introduce canonically normalized fields — (24) 1//2 $ 
and r^v = (24)~ 1 / 2 r. Then the Lagrangian (76) becomes 

^2,scalar[ ( i?Af(0: Tjv(^)] = 

]4> 2 N + f jv$jv - ^m%{^% + 2T N <P N + aT 2 N ) , 



where as usual ttiq = 61m 2 , and 

= 3(bi+26 2 ) 
01 ~ 2(b 1 + b 2 ) 

Going to Fourier space, we get 



S- 



1 



2, scalar 



Z7T 



,f_ a) )M(o;) 



with 



M(oj) = 



itIq —arriQ 



(A8) 



(A9) 



(A10) 



(All) 



Inverting this quadratic form we find that the propagator 
has poles at w 2 = 



and at 



-u> 



(h + 46 2 ) 

i 

2(6l+62) 



(A12) 



which therefore are the masses of the two degrees of free- 
dom in the scalar sector. The negative sign in front of u> 2 
in eq. (A12) reflects the fact that this mode is ghost-like, 
hence its kinetic term has the opposite sign compared to 
a standard kinetic term. 
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Appendix B: Derivation of j° — 0, dif = from the 
non-dynamical equations of motion 

As we mentioned in Sect. IV, defining j" as in eq. (131), 
the use of the non-dynamical equations of motion ensures 
that j° = and dij 1 = 0. To prove this result we write 
j v using the harmonic decomposition of the metric. For 
the temporal component we get 



3 



-(V 2 7 + 60). 



(Bl) 



Using cqs. (84) and (86) we find j° = 0. Therefore, the 
condition (67) with v = is indeed a consequence of the 
non-dynamical equations of motion. 

For the spatial component j % the situation is different. 
In harmonic variables, 



-0*7 - 20 4 + \d l V 2 \ + \v 2 e l - ft 



-26 2 (0*V> + 30*0) . 



(B2) 



This quantity depends both on the scalar functions 
-0,0, 7, A, and on the vectors e 4 and /3 4 . The part de- 
pending on the scalar functions vanishes upon use of the 
non-dynamical equations of motion (84) and T = 0. We 
remain with 



V 



2/3' 



(B3) 



We found in eq. (107) that indeed the equations of motion 
in the vector sector imply the condition V 2 e* — 2/3* = 0. 
However, to get this result we used the dynamical equa- 
tion of motion (106) obtained performing the variation 
with respect to e\ Thus, the use of the non-dynamical 
equations of motion is not sufficient to derive the condi- 
tion = 0. For this, it is necessary to use the full equa- 
tions of motion, including a dynamical component. How- 
ever, e 1 and f3 l are transverse vectors, so from eq. (B3) 
it follows dij 1 = 0. This means that, at the FP point, 
the non-dynamical equations are sufficient to derive the 
conditions 



f = o, o,y o. 

and therefore also d^j^ = 0. 



(B4) 



respect to ip,j,\ and <j> one finds 
1 



Sip : (d-1)- 



■V 2 $ 



d(d- 1)H($ + HV 



#7 



2(d-l) 



7, 



(d-2)$-* = ^A ; 



<I> 



dH<£ 

2 



H$ + dH 2 ^ 



1 



[(d-i)<f> + i/>]. 



-dm 2 g (j) , 
(CI) 

(C2) 

(C3) 
(d-2)*] 
(C4) 



Using cqs. (C2) and (C3) one can eliminate 7 and A, and 
then the definitions of $ and eqs. (177) and (178), al- 
low us to express 0, ip in terms of and \I r . Plugging the 
resulting expressions into eqs. (CI) and (C4) and combin- 
ing the equations one finds the coupled set of equations 



1 



V 2 ($ + #) + dH($ + *) + d(d - 3)H 2 ^ 

+ [2d(d - 2)H 2 - dm 2 ] $ = , (C5) 



and 



[n-m 2 g ) [a- 2 (<S> + *)] =0, 



(C6) 



where □ = (1/y/— g)d^(y/— g g^d u ) is the dAlembertian 
in curved space on scalar functions. 14 When specialized 
to d = 3, cqs. (C5) and (C6) agree with the result pre- 
sented in [27], after the appropriate change of notation. 15 
In flat space eqs. (C5) and (C6) reduce to 



V 2 ($ + *) = dm 2 g $>, 



($ + *) =0. 



(C7) 
(C8) 



Applying the Laplacian to eq. (C8) and using eq. (C7) 
we get 



m 2 (D 



(C9) 



and therefore, for m g ^ 0, we recover our flat-space result 
(83). 

At first sight eq. (C6) seems to imply that, in de Sit- 
ter, the radiative degree of freedom in the scalar sector is 
described by a _2 ($ + The situation is however more 
subtle. Equations (CI) and (C4) have been obtained by 
treating {0, -0,7, A} as independent fields. At this level 



Appendix C: Equations of motion in terms of $ and 
* in de Sitter 

In this appendix we compare our analysis with the dis- 
cussion of the equations of motion in de Sitter performed 
in [27]. Taking the variation of the action (182) with 



In the derivation one uses the fact that in de Sitter, for any scalar 
function /, we have -a 2 n(a" 2 /) = f + (d-A)H f -2(d-2)H 2 f - 
a- 2 V 2 /. 

In particular, in [27] our <3> is called and our \t is called <3>, the 
metric signature is the opposite and the equations are written in 
conformal time. See also [21] for similar perturbation equations 
in the context of bigravity. 
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$ and \& are therefore simply a notation for the combina- 
tions given by eqs. (177) and (178). In particular, since 
^ includes a terms A, eq. (C3) is a fully dynamical equa- 
tion for the variable A, and not a constraint that allows 
us to eliminate A from the theory in favor of ^> and $. A 
way to understand this point is to observe that, in order 
to have a well-defined Cauchy problem, we must assign 
the values of hij and at some initial time to. In the 
scalar sector, this means that we assign {</>(i ), A(£o)} and 
{<^(io)> A(io)}- However, for H ^ 0, $ contains a term 
A and VP contains terms A and A. Therefore, assigning 
A (to) and A(io) is not sufficient to provide the initial val- 
ues ($ + #)(t ) nor ($+*)(£ ). Thus, eqs. (C7) and (C8) 



arc not a closed set of equations that can be solved for 
$ and ^P, once given the initial conditions on the metric 
(in which case one could have then determined A from 
eq. (C3)). We still need eq. (C3), written as a second- 
order differential equation in A, to evolve the system. 

Thus, these manipulations of the equations of motion 
are not the correct way of eliminating the non-dynamical 
variables, and ($ + is not the radiative degree of free- 
dom in the scalar sector. In order to correctly identify 
the field that describes the radiative degree of freedom 
in the scalar sector we must go through the procedure 
presented in sect. VII. 
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